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The two-dimensional (2D) Heisenberg magnet with a weak 
easy-plane anisotropy is considered. A renormalization group 
(RG) analysis in this model is performed for both quantum 
and classical cases. A crossover from the Heisenberg to 2D 
XY model is discussed. The magnetic transition owing to 
the interlayer coupling is considered. Analytical results for 
the Kosterlitz-Thouless and Curie (Neel) temperatures are 
derived with account of two-loop corrections. The results are 
compared with experimental data, e.g. on K2CUF4, and turn 
out to provide a quantitative description, unlike the standard 
one-loop results. 



The interest in the magnetic properties of layered sys- 
tems has been recently greatly revived. It is well known 
that even weak magnetic anisotropy can play in such sys- 
tems an important role. In the present paper we discuss 
the case of the easy-plane localized-spin system. The sim- 
plest classical two-dimensional (2D) XY model was stud- 
ied in detail [HQ, and the relevance of topological (vor- 
tex) excitations in thermodynamics was established. In 
particular, the Kosterlitz-Thouless transition, connected 
with unbinding of the vortex-antivortex pairs, was found. 
The transition temperature, where power-law behavior of 
spin correlation function is changed by exponential one, 
is estimated as 



-JS 2 



(1) 



where J > is the exchange integral. In the quan- 
tum XY model the situation is still more complicated, 
since not only transverse, but also z-components of spins 
should be taken into account. 

A different situation takes place in both quantum and 
classical Heisenberg 2D model with a weak easy-plane 
anisotropy, which is a more physically real case Q3|. A 
simple expression for Kosterlitz-Thouless temperature 
obtained in Ref. B reads 



T, 
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KT 



ln[7rV(l - V)] 
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(r) = J z /J x - y is the anisotropy parameter) and has the 
same form as the result for the magnetic ordering point 
of the easy-axis layered magnet 01. As well as the lat- 
ter result (see discussion in Refs. g-Q]), the formula (||) 
is insufficient for a quantitative description of the ex- 



narrow region near Tkt- The situation is reminiscent 
of the easy-axis layered magnet where the topo- 

logical (domain-wall) excitations are important only in 
the vicinity of the Curie (Neel) temperature Tc(Tn) In 
such a situation, similar to Jt| , the renormalization-group 
(RG) analysis can be performed to calculate Tkt with 
higher accuracy. This analysis is the aim of the present 
paper. Further on, we consider effects of interlayer cou- 
pling, which lead to occurrence of the true long-range 
magnetic ordering, and calculate Tq(Tn). 

We consider the easy-plane Heisenberg model 



[SfS? 
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where Jy = J ( J > in the ferromagnetic (FM) case and 
J < in the antiferromagnetic (AF) case) for the nearest- 
neighbor sites i,j in the same plane, and Jij = a J for 
i,j in different planes, 77 < 1. We suppose 1 — 77, a <C I . 
Note that the effect of the single-site anisotropy 



h d = dJ2(sD 2 , d>o 



(4) 



is the same as that of exchange anisotropy with I — rj = 
D(l - 1/25) /4|J| provided that D/ \J\ < 1. 

The partition function of the model (||) can be repre- 
sented in terms of a path integral over coherent states 
(see, e.g., Refs. §,§): 

Z = J Dtt5{tt 2 - 1) exp(-£ dy „ - C st ) (5) 

1/T 

£dyn = iS 1 J A(7Tj)-^ 



1/T 

Art = ^S 2 j (IT Jij [^Xi^Xj + KyiTTy 



Tyj + 1Wzi^zi\ 



where 7r is the unit-length vector, and A(tz) is the vector 
potential of the unit magnetic monopole, which satisfies 
the equation V x A(tt) • 7r = 1. 

In the classical case (i.e., with £dyn being neglected) 
we have two types of excitations: the field tt v describes 
the gapless in-plane excitations, and the field ir z describes 
the out-of-plane excitations with a gap. Expanding (|5|) 
in 'Ky.z^ x being a (local in AF case) spin quantization 



perimental data (see Ref. ||). Since T KT « JS 2 , one axigj ' Q leading Qrder in y s we haye 



can expect that thermodynamic properties of these sys- 
tems are determined by usual spin waves, except for a 
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where Q is the wavevector of magnetic structure. The dy- 
namical part Cdyn, which is present in the quantum case, 
results in (i) quantum renormalizations of the Hamilto- 
nian parameters (in AF case only), which are supposed 
to be already performed and (ii) the summations over 
wavevectors are bounded by \JT j JS in FM case or T/c 
in AF case (c is the quantum-renormalized spin- wave ve- 
locity) rather than by the Brillouin zone boundary (see 
Ref. §). 

The interaction of spin waves, which occurs in higher 
orders in l/S 1 , leads to temperature renormalizations 
of the Hamiltonian parameters. Due to the smallness 
of anisotropy, large logarithms occur in these renor- 
malizations: ln|T/(l — rj)J\ in the quantum FM case, 
ln[T 2 /(l— 77) J 2 } in the quantum AF case, and ln[l/(l— 77)] 
in the classical case. It is natural to sum up these log- 
arithms within the RG approach. However, there ex- 
ists the difficulty owing to that gapless ^-excitations 
are present too. In the absence of the interlayer cou- 
pling, they lead to infrared divergences in some quan- 
tities like the in-plane magnetization. In the presence 
of interlayer coupling, another type of logarithms occur, 
ln(T/aJ), ln(T 2 /aJ 2 ) or ln(l/a), depending on the case. 
The situation, where two types of excitations with differ- 
ent characteristic scales are present, is typical for systems 
demonstrating a crossover |lC| , In our model this is the 
crossover from the Heisenberg (almost isotropic) behav- 
ior to the XY behavior. 

To describe correctly this crossover we include 
anisotropy in all the renormalization factors p0| . Wc 
also introduce the scaling factors of the field 7r. Be- 
cause of anisotropic character of the model, we have 
two such factors: Z xy and Z z , so that tTxr/t^x = 
^vrI^v — Zxy and tTzr/kz — Z z . We use the nor- 
malization condition T zz (0) = 1 — 77 (which fixes 
the gap of z-excitations) instead of the standard one, 
dT z 2 }{q)/d{q 2 ) — 1, T Z 2 J (q) being the two-point vertex 
function (inverse Green's function) of the field tt z . Then 
we have Z z = 1. For other Hamiltonian parameters we 
obtain the following system of RG equations 
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(7c) 
(7d) 



M^ = 2t„/fo M ,M)+0(£) 



d\noe u , 



where /j, is the scale parameter, /(?7 M , /i) = r/^fj, / (r/^^i 2 

i-v), 



T/(2ttJS 2 ) FM 
T/(2tt Ps ) AF 



is the dimensionless temperature, p s ~ S(S + 0.079) |J| 
being the spin stiffness pj] |. First two equations in (Q) 
are written down to two-loop order, while last two to one- 
loop order, which is sufficient to obtain the final results 
to the two-loop order accuracy. The initial scale /z for 
these equations is 



^0 



'32 
T/c 



classical regime (T ^> \J\S) 
quantum regime (T <C I^IS*), AF , 



y/T/JS quantum regime (T < JS), FM 



for the details see, e.g., Ref. 0. 

The flow of RG parameters is shown schematically in 
Fig.l (for comparison, the easy-axis case 77 > 1 is de- 
picted too). From the equations (7b) and (|7q) we obtain 
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where the function <&(/!) = 0(t^) comes from 0-terms in 
(0) and corresponds to the contribution of higher-order 
loops. For fi ^> y/1 — 77 the effective temperature is 
small (which guarantees that the spin- wave theory works 
well), so that we have $(/z) <§C 1 and 



For fi <C — 77 we obtain 
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(10) 



and in this regime is //-dependent only through the 
function <&(//). The scale l/^l — 77 is just a character- 
istic scale for the crossover from the Heisenberg to XY 
behavior and (|hJ) describes < M in the XY regime. On 
the other hand, in this regime only vortices contribute to 
the temperature renormalization since such a renormal- 
ization owing to spin-waves is absent (the interaction of 
spin waves in the XY model is due to topological effects 
only). Thus for the temperature renormalization we have 
the system of RG equations |l|,0| , which in our notations 
can be written as 



dy, 



(11a) 
(lib) 



It should be noted that the coupling constant for the 
vortex system is not t (as for spin waves), but y = 
exp(—Eo/T) where Eq is the energy of a vortex core. 
Therefore equations (O) are applicable for small enough 



fi. Let /ii <C \/l — 77 be the scale where we pass to scaling 
(|TT|). Then the solution of Eqs. ( pi] ) for t > tur reads 
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where 



= 4 + 2Ci tan ( C x In — + C* 2 
Mi 



Ci - V(8^/i*i) 2 - (4*i - l)7(2ii) 
1 - 4ti 



(12) 



tanC*2 = 



^(STrj/iii) 2 - (4*i - l) 2 



and *i = * M1 , ?/i = y M1 arc determined by 



*i * vT^J? 
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(14) 



It should be stressed that even if the original model is 
quantum one, the resulting XY model is classical since 
Hi < V^-V < -^r 1 ( L r = JS/T for FM case and L T = 
c/T for AF case is a characteristic length for quantum 
effects) and at scales much larger than L T quantum and 
classical systems becomes indistinguishable. Thus all the 
quantum effects are already taken into account at the 
scales /i 3> \/l — rj, where the behavior of RG trajectories 
is Heisenberg one. 

The Kosterlitz-Thouless temperature Tkt is deter- 
mined by the equation of the separatrix line for Eqs.(|T] 



87ryi = l/*i - 4, * = t KT 



(15) 



This line separates the low- and high-temperature phases. 
For small enough /i we have $(^) — >const, d$(fi)/dfi — > 
and we have for tKT = Tkt/(2ttJS 2 ) (or Tkt/(2kPs) in 
AF-case) 



tKT 



ln^o/VT - Tj) + 2 1h(2/* K t) + C 



(16) 



where C = 4 — 6 In 2 — $(/i — > 0) is an universal constant. 
This result is identical with that for the Curie (Neel) 
temperature of an easy-axis magnet @, except for the 
constant C, which needs not be the same as for the easy- 
axis case. 

In the critical region above tKT, 



i-(*-i T -*-i)<l, 



(17) 



the expression for the correlation length obtained from 
(O) reads 



Mi yfi - V 



A 



exp 



2\lfg T -t- 1 



(18) 



and has the same form as for the XY model (A is a 
constant). Under the condition, opposite to (17), we have 
the standard Heisenberg behavior jll| 



£ = (Q//i )*exp(l/t) 



(19) 



In the presence of interlayer coupling, the magnetic or- 
dering at low enough temperatures occurs. Due to topo- 
logical effects, the transition temperature grows up from 
Tkt% and not from zero. Note that in this case Tkt 
plays a role of a crossover temperature from 2D to 3-D 
XY behavior rather than a critical temperature, and the 
only true phase transition is connected with the magnetic 
ordering at Tc(Tjv). 

In the case a <C 1 — 77 we choose /_ti such that 
a 1 / 2 < |Ji < (1 - v) 1 ^ 2 - In terms of RG transforma- 
tion (see Fig.l), at fi = /ii we have not 2D, but quasi-2D 
XY effective model with the lattice constant Mo/Ml an d 

'0 



the interlayer coupling (/jq/mi) 2q: i where, as follows from 



Oil 



at/ti 



(20) 



With further flow of the RG transformation we should ar- 
rive at the 3D XY model. However, this part of the RG 
transformation meets with difficulties owing to a compli- 
cated geometry of vortex loops (see Ref. [[L2) and refer- 
ences therein). Instead of direct calculation of RG tra- 
jectories, we use the same scaling arguments as in Ref. 
Q. The transition temperature can be estimated from 
the requirement that the correlation length of the model 
without interlayer coupling (a — 0) coincides with the 

characteristic scale of the crossover from 2D to 3D XY 

1/2 

model, 1/a/ (in the units of the lattice constant of orig- 
inal lattice). Then we have for the critical temperature 
T c /(2nJS 2 ) (or T N /{2np s )) in the case a < 1 - 77 



tr. 



Ill ■ 



Mo 



2 hi- 



tKT 



C 



A 2 



(21) 



The last term in the denominator determines the differ- 
ence between * c and *xt- Since this term can be not too 



small, we do not expand (21) in it. 

The result (^lj) is qualitatively valid up to a of order of 
1 — 7] (in this case last term in the denominator leads to 
renormalization of C only) . Consider now briefly the case 
a ^> 1 — 77. Then the corrections to the RG result for the 
quasi-2D magnets owing to the easy-plane anisotropy 
are given by 



tr. = 
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2 
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(22) 



where %j) — ^3(2 — 7^) is the crossover exponent, 1/3 is 
the corresponding critical exponent for the 3D Heisen- 
berg model and 7^ is the anomalous dimensionality of 
the anisotropy parameter near 3D Heisenberg fixed point, 
see, e.g., Ref. [ |l0fl . The e-expansion in the anisotropic 
4 — e dimensional </> 4 model for e = 1 (which has the 
same symmetry as the model under consideration) yields 
ip ~ 0.83, see Ref. ]K(| . For an antiferromagnet, the con- 
stant C ~ —0.066 was calculated within the 1/N expan- 
sion ||. Unlike (j2l|), the last term in the denominator of 



3 



( p2| ) has not inverse-logarithmic form. This is a conse- 
quence of the fact that the correlation length in the 3D 
Heisenberg model does not demonstrate the exponential 
behavior (1/3 is finite) . By this reason the correction in 



the denominator of (22) is small and can be neglected. 

Finally, we consider the experimental situation for lay- 
ered magnets. The mostly investigated easy-plane sys- 
tem is the compound K2CUF4. This is a S — 1/2 ferro- 
magnet with Tkt = 5.5K, Tc = 6.25K and the param- 
eters J = 20K, 1 - 7} = 0.04, a = 6 • 10" 4 (see, e.g., 
Ref. ||). Substituting these values into (|l^) and ( pil| ) we 
obtain C ~ —0.5 and A ~ 3.5. Note that the formula (||) 
yields the value Tkt = H-4K which is much larger than 
the experimental one. 

Another example of a quasi-2D FM AY"-like system 
is the stage-2 MCI2 graphite interlayer compound with 
5 = 1. According to Ref. §j, J = 20K, 1 - 7? = 8 • 
10~ 3 and a — 5 • 10~ 5 . Using the same values of A 
and C as for K 2 CuF4, we calculate Tkt = 17.4 K and 
Tc = 18.7 K, which is in agreement with experimental 
data (both values Tkt an d Tq lie in the region 18 — 20 
K). At the same time, using the formula (||) yields Tkt = 
35. 3K, which is again twice larger as compared to the 
experimental value. 

We have also applied our results to the compound 
BaNi2(PC>4)2 which is a S — 1 antiferromagnet with 
I J I = 22.0 K and easy-plane anisotropy 1— 77 = 0.05, a — 
1 • 10~ 4 , see Ref. [§. We obtain T KT = 23.0 K which 
coincides with the experimental value and Tjv = 24.3 
K, again in excellent agreement with T^ xp = 24.5 ± IK. 
Note that in spite of Tkt ~ \J\S for this compound, the 
true criterium of the quantum regime is (T/JS) 2 <C 32 
(see Ref. 0), and this case also should be considered as 
quantum one. 

To conclude, we have investigated the Heisenberg 
model with a weak easy-axis anisotropy. We have per- 
formed the two-loop RG transformation with unknown 
function (which takes into account the contribu- 

tion of higher loops and non-spin-wave excitations) and 
joined the results with well-known behavior of the RG 
trajectories in the 2D XY model. In such a way we have 
obtained simple analytical expressions for the Kosterlitz- 
Thouless and Curie (Neel) temperatures. These expres- 
sions contain two constants A, C which are still inde- 
terminate within the RG approach. The calculation of 
these constants, as well as of the corresponding parame- 
ters for the isotropic quasi-2D and easy-axis 2D Heisen- 
berg model ||], is possible by numerical (e.g., by the 
quantum Monte-Carlo) methods |u|. At the same time, 
our results already enable one to estimate the Kosterlitz- 
Thouless temperature (and also to determine the differ- 
ence Tc(Tjv) — Tkt) with the accuracy which is sufficient 
to fit experimental data on layered magnets, unlike the 
simplest expression (||). 

We are grateful to B.N.Shalaev for useful discussions. 

Figure caption 

Schematic picture of the RG trajectories in layered 
magnets. Left-hand side: the flow from the 2D easy-axis 



Heisenberg (H+EA) to 2D Ising model. Right-hand side: 
the flow from the 2D easy-plane Heisenberg (H+EP) to 
2D XY model. The inflection points c\, C2 mark the 
crossover regions. The dashed lines are for the corre- 
sponding quasi-2D models. 
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